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The measurements of the mean lifetime of the r lepton are reviewed. The conditions for measuring the lifetime 
at various e + e~ colliders are compared and the analysis methods are briefly described. The new developments 
since the previous Workshop on Tau Lepton Physics are listed. The world average is t t — 290.5 ± 1.0 fs. The 
LEP experiments dominate this average and have analyzed nearly all of their data. In anticipation of the next 
era of precision measurements at CLEO and the b factories, the important sources of systematic errors and the 
treatment of systematic biases are discussed. 
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1. Introduction 

The mean lifetime of the r lepton is deduced 
from geometrical reconstruction of r daughter 



tracks in e + e 



t+t events. We need to 



know the r lifetime in order to make certain tests 
of lepton universality. Such tests are sensitive 
to the tj/ t W coupling and also to possible new 
physics 0]; at present the precision is limited by 
the experimental uncertainties on the r lifetime 
and branching fractions. The r lifetime is also 
useful for evaluating the strong coupling constant 

as I- 

In this talk I use the conventional coordinate 
system for e + e _ experiments, i.e., with the z 
(polar) axis along the direction of the incident 
beams. The impact parameter d of a recon- 
structed charged track is measured in the xy 
projection with respect to the nominal interac- 
tion point; d is signed according to the z com- 
ponent of the track's angular momentum about 
this point. The azimuthal decay angle ip of a r 
daughter track is defined as the signed quantity 
^daughter ~ 4>t in the laboratory frame, where ip 
lies between — 7r and 7r. 

At present, useful measurements of the r life- 
time can be obtained by SLD, the four LEP ex- 
periments, and CLEO, and the experimental con- 
ditions are quite different in these places. As far 
as the lifetime measurement is concerned, the rel- 
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evant differences include the r flight distance in 
the laboratory (which varies as f3j of the r boost) , 
the opening angles of the r decay products (which 
vary roughly as 1/7) , the dimensions of the lumi- 
nous region, and the number of collected r pairs. 
On the other hand, the impact parameters of the 
r daughters are on the order of ct t = 87 /im in all 
experiments, as long as the r has /3 ~ 1. Com- 
parisons of the experimental conditions are given 
in Table | and Fig. |. 

2. Methods 

In this section I give a brief description of the r 
lifetime measurement methods that are currently 
in use. 

2.1. Decay length method 

The decay length (DL) method (or vertex 
method) is applied to three-prong r decays ||. 
The procedure consists of reconstructing the r 
decay vertex in two or three dimensions. In or- 
der to evaluate the r decay length, an estimate 
of the r production point is also needed. In an 
e + e~ — ► t + t~ event, the only available estimate 
is the centroid of the luminous region. Because 
the luminous region is huge along the direction of 
the beam axis (considerably larger than the typ- 
ical r flight distance), we can effectively estimate 
only two of the three coordinates of the r produc- 
tion point, and we end up with a measurement 
of the r displacement L xy in the xy projection. 
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Table 1 

Experimental conditions for t lifetime measure- 
ments. Here, N TT is the approximate number of 
produced e + e~ — > t + t~ events in the data sam- 
ple, E T is the t energy in the laboratory (neglect- 
ing radiative effects), 7 refers to the boost of the 
r, and (3jct t is the mean r flight distance in the 
lab. 





SLD 


LEP 


CLEO 


Nrr 


20K 


4 x 200K 


5M a 


E T (GeV) 


45.6 


45.6 


5.3 


7 


25.7 


25.7 


3.0 


(3jct t (jum) 


2200 


2200 


240 



a Collected with the present (CLEO 2.5) detector 
configuration. 



An estimate of the polar angle 8 T of the r di- 
rection, taken for example from the event thrust 
axis, is therefore needed in order to calculate the 
t flight distance L in three dimensions. The mean 
lifetime is deduced from the mean of the L dis- 
tribution and the mean value of ffy; the latter is 
calculated from simulated t + t~ events, including 
initial and final state radiation, after the event se- 
lection criteria are applied. 

Modern vertex detectors can measure precise 
impact parameters in both the r<f> and rz views. 
Because the luminous region is so large along z, 
the possibility of measuring the z coordinate of 
the r decay vertex is not particularly useful for 
the classical DL method. However, it is impor- 
tant to realize that the tracking information from 
the rz view can also significantly improve the 
measurement of the (x, y) coordinates of the de- 
cay vertex, from which the lifetime is extracted. 
This point is illustrated in Fig. ||[ 

At LEP and SLD, the size of the luminous re- 
gion and the tracking resolution are such that the 
statistical uncertainty on the mean decay length 
is dominated by the natural width of the exponen- 
tial t distribution; the relative uncertainty on the 
mean lifetime is not far from its optimum value, 
Ar T /r T = 100%/v^Vrj where N T is the number 
of selected t decays. The systematic errors in DL 
analyses at high energies also tend to be fairly 
small. In short, our detectors and our technique 
for measuring the r lifetime from three-prong de- 
cays at LEP and SLD are very effective. 

At lower center-of-mass energies the r decay 
length is shorter and the size of the luminous re- 




Figure 1. Experimental conditions at SLD, LEP, 
and CLEO. All drawings show the same three- 
prong r decay, boosted to the appropriate energy 
for each experiment and projected onto the xy 
plane (20 x actual size). The ellipses represent- 
ing the luminous region include the typical un- 
certainty on the beam axis coordinates. 



gion can significantly dilute the precision of the 
measurement. In such cases the decay length res- 
olution can be improved by considering the sep- 
aration of the two decay vertices in 3-3 topology 
events. This procedure eliminates the smearing 
due to the size of the luminous region, but it does 
not dramatically improve the final results because 
the statistics are low and the qq background is 
larger in the 3-3 channel. 

2.2. One-prong decays 

While the DL method is entirely satisfactory 
for analyzing three-prong r decays, we cannot ap- 
ply such a straightforward technique to the one- 
prong decays. Due to the unobserved neutrinos, 
the t direction is unknown. We cannot recon- 
struct the decay length of an individual r decay- 
ing into one prong. We can nevertheless measure 
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Figure 2. At decay in which the three charged 
daughter tracks emerge in the same azimuthal di- 
rection ip. The vertex cannot be reconstructed 
from the track measurements in the r<f) view 
alone, but a full three-dimensional vertex recon- 
struction is possible when the measurements in 
the rz view are added. 



the mean decay length from a collection of one- 
prong decays. The relative statistical uncertainty 
will be considerably larger than 100%/y/N T . We 
will do best if we analyze both r decays in an 
event together. 

It is difficult to incorporate all of the event in- 
formation into a simple method. We now use sev- 
eral methods to analyze one-prong decays. Each 
method uses a different subset of the available in- 
formation, and none of the methods is vastly su- 
perior to the others. We combine the results from 
the various methods, taking into account the sta- 
tistical and systematic correlations, to utilize as 
much of the available information as possible. I 
now proceed to describe these methods in more 
detail. 

2.3. Impact parameter method 

The impact parameter (IP) method is applied 
to one-prong r decays [0. In this method, an 
estimate of the r direction is taken, for example, 
from the event thrust axis. The lifetime-signed 
impact parameter of the daughter track is then 
defined: 

n= J d if^>0; 
\ -d if V < 0. 

The mean of the D distribution is then roughly 
proportional to r r . The dependence of the D dis- 
tribution on t t is determined from Monte Carlo 
simulation. 

Decays with D < result from the impact pa- 
rameter resolution, the size of the luminous re- 
gion, and errors on the r direction. The last of 



these effects is probably the most dangerous be- 
cause it brings about a substantial change in the 
mean of D, and we rely on the Monte Carlo to cor- 
rectly describe the r direction errors. Although 
we tend to think of the IP method as a one-r-at- 
a-timc method, the r in the opposite hemisphere 
contributes to the thrust axis determination and 
hence affects the D distribution. 

2.4. Impact parameter difference method 

The impact parameter difference (IPD) method 
is applied to 1-1 topology events ||. In this 
method the mean r decay length is extracted by 
considering the correlation between the difference 
on the daughter track impact parameters d and 
the difference of their azimuthal angles cf>. Specif- 
ically, we define Y = d+ — d- and X = A<fi sin 9 T , 
where Acf> = 4> + — <f>- ± 7r is the acoplanarity of 
the two daughter tracks. If the r + and t~ are 
back to back in the xy projection and the decay 
angles %j) are small, such that sin^ = tp, we find, 
at a particular value of A, that (Y) = LX, where 
L is the mean r decay length in the lab. A fit to 
the Y vs X distribution is performed to extract 
the slope L. The polar angle 8 T is taken from the 
event thrust axis; the resulting error on the r di- 
rection has a negligible effect on the fitted L. The 
main disadvantage of the IPD method is that the 
uncertainty on the t + t~ production point due to 
the size of the luminous region enters twice in the 
smearing on Y. 

2.5. Impact parameter sum methods 

The miss distance (MD) U and momentum- 
dependent impact parameter sum (MIPS) @ 
methods are designed to give improved statisti- 
cal precision by virtually eliminating the smear- 
ing effects related to the size of the luminous re- 
gion. In a 1-1 topology event we define the "miss 
distance" A = d + + d- . This sum of impact pa- 
rameters is, roughly speaking, the distance in the 
xy projection between the two daughter tracks 
at their closest approach to the beam axis. This 
quantity is almost independent of the r + r _ pro- 
duction point. The A distribution depends on r r ; 
a Monte Carlo simulation is used to parametrize 
the true distribution, in order to extract the life- 
time from the data. The main disadvantage of 
these methods is that the results of the fit to the 
data are sensitive to the assumed impact param- 
eter resolution. 

I refer to the simplest form of this analysis as 
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the MD method. The MIPS method is a re- 
finement of MD in which the A distribution is 
parametrized in terms of the momenta in the 
lab of the two r daughter tracks. DELPHI'S 
new results announced at this workshop feature 
the MIPS refinement and one other: the A dis- 
tribution is parametrized separately for lepton 
and hadron daughters; I refer to this method as 
MD++. 

2.6. Three-dimensional impact parameter 
method 

The three-dimensional impact parameter 
(3DIP) method makes use of more of the kine- 
matic information in the events to yield a higher 
sensitivity per event than the other one-prong 
methods (|]. The main disadvantage is that the 
method can only be applied to hadron vs hadron 
events (42% of all r+r~ events). Because a 
t — > hadron decay yields only one unobserved 
neutrino, it is possible to reconstruct the r direc- 
tion in hadron vs hadron events up to a twofold 
ambiguity. Let fi and T2 denote the two possi- 
ble t" directions reconstructed for a particular 
event (Fig. ||). If we then project the event along 
a direction chosen such that T\ and T2 coincide, 
we end up with no uncertainty on the r direction 
in that projection. We then define a generalized 
impact parameter sum in that projection, so that 
there is almost no smearing due to the size of 
the luminous region. A fitting procedure oper- 
ates on this impact parameter sum and on the 
two projected r decays angles in order to extract 




Figure 3. In the 3DIP method, the event is pro- 
jected along the direction given by fa — fi, where 
fi,2 are the two possible t~ momentum direc- 
tions. 



the mean lifetime. 

The 3DIP method is the first to use impact 
parameter information from the rz view in the 
analysis of one-prong decays. The method has the 
extremely important advantage that the tracking 
resolution and the r lifetime can be extracted si- 
multaneously from the t + t~ events. 

3. New lifetime results since TAU96 

There are four new developments to report: 

• At TAU98, L3 is reporting preliminary re- 
sults from their 1994-95 data, analyzed with 
the DL and IP methods [§. The new L3 
average (1991-95 data) is r T = (291.7 ± 
2.0[stat]±1.8[syst])fs. 

• At TAU98, DELPHI is reporting preliminary 
results from their 1994-95 data, analyzed 
with the DL, IPD, and MD++ methods @. 
The new DELPHI average (1991-95 data) is 
t t = (291.9 ± 1.6 [stat] ± 1.1 [syst]) fs. 

• The thesis of Patrick Saull (ARGUS) f| de- 
scribes the vertex impact parameter (VIP) 
method, which provides improved lifetime 
sensitivity for 1-3 topology events in cases 
where the size of the luminous region limits 
the precision of the DL method. The VIP 
method uses the impact parameter of the 
one-prong track with respect to the three- 
prong vertex, and the acoplanarity of the 
one- and three-prong jets. (A similar ap- 
proach is described in Jl2[ .) 

• In 1997, ALEPH published results from the 
3DIP method (1992-94 data) § and from 
the MIPS, IPD, and DL methods (1994 
data) [p"3f , preliminary versions of which had 
been shown at TAU96. 

4. Summary of measurements 

In calculating the world average r lifetime, I 
follow the Particle Data Group |Q and ignore 
early measurements with large uncertainties. The 
measurements are listed in Table || and plotted in 
Fig. [|. In most cases the results shown are them- 
selves averages of two or more measurements ob- 
tained by a given experiment with different meth- 
ods and/or data samples. 

The world average is r T — 290.5 ±1.0 fs, where 
the systematic errors in the various experiments 
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Table 2 

Measurements of the r lifetime. 



Experiment 



ALEPH 
DELPHI 
L3 @* 
OPAL |X5| 



PT 



t t ± stat ± syst (fs) 



CLEO II |Q 



SLD |17| 




290.1±1.5±1.1 
291.9±1.6±1.1 
291.7±2.0±1.8 
289.2±1.7±1.2 

289.0 ±2.8 ±4.0 

288.1 ±6.1 ±3.3 



* Preliminary 



are assumed to be uncorrelated. The x 2 describ- 
ing the consistency of the measurements is 1.36 
for 5 degrees of freedom, corresponding to a con- 
fidence level of 0.929. The four LEP experiments 
contribute 94% of the total weight in the average. 
Since the beginning of the LEP era, the uncer- 
tainty on the world average has been reduced by 
a factor of 8. By now, almost all of the LEP1 
data has been analyzed. (The LEP2 data is not 
expected to yield any useful r lifetime results.) 

5. Systematic errors 

Most recent measurements of the r lifetime are 
statistics limited. Nevertheless, it is useful to ex- 
amine the systematic effects that we will need to 
deal with in the coming years in order to further 
improve the precision of the measurements. Some 
of the important sources of systematic errors are 
tracking errors (simulation and/or parametriza- 
tion); vertex fit and lifetime extraction proce- 
dures; and detector alignment. I will now discuss 
each of these topics in turn. 

5.1. Tracking errors 

I would like to mention two delicate issues re- 
lated to tracking errors. 

The first issue concerns the dependence, in 
some methods, of the measured lifetime on the 
assumed impact parameter resolution. In such 
cases, it is mandatory to measure the resolu- 
tion from reconstructed tracks in the real data. 
Bhabha and dimuon events are readily available 
for this job, but the high momentum electrons 
and muons in those samples are not representa- 
tive of the r daughter tracks and the contribution 
to the impact parameter resolution from multiple 
scattering cannot be studied. Some experiments 
employ 77 — > e + e~, events to parametrize 

the resolution at the low end of the momentum 



range. While these test samples can give a fairly 
precise description of the impact parameter res- 
olution for electrons and muons, it is not easy 
to use the real data to parametrize the effects of 
nuclear interactions on pion and kaon tracks. 

The second issue concerns the correlation be- 
tween the errors on the reconstructed impact pa- 
rameter and direction of a track, e.g., between 
d and <f>. This correlation is positive and results 
from the extrapolation of the reconstructed tracks 
from the measured points in the tracking detec- 
tors to the interaction region. In some meth- 
ods the correlation can simulate a longer r life- 
time jlj|. The effect is especially bad for the 
t (compared to other particles) because (1) the 
short r lifetime leads to small impact parameters, 

(2) the small t mass leads to small decay open- 
ing angles (which get smaller at higher t/s), an d 

(3) we fit to the entire proper time spectrum (as 
opposed to the situation in charm lifetime mea- 
surements in fixed target experiments, where a 
cut L > L m i n is imposed and the mean lifetime 
is determined from the slope of the proper time 
distribution). The effects of the tracking errors 
on the measured r T must be carefully taken into 
account. 

5.2. Vertex fit and lifetime extraction 
procedures 

Although the direct reconstruction of r decay 
lengths in the DL method appears to be quite 
straightforward, several subtle effects are present, 
yielding biases on the measured lifetime. These 
effects are related to the tracking resolution, and 
they can be substantial in experiments where the 
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Figure 4. Measurements of the r lifetime. 
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vertex resolution along the r direction is larger 
than the mean decay length (not the case at SLD 
and LEP). 

• Due to the correlation between the track im- 
pact parameter and direction errors, fluctu- 
ations to larger opening angles in a three- 
prong decay tend to be associated with up- 
ward fluctuations in the reconstructed de- 
cay length. The larger opening angles also 
lead to a smaller calculated uncertainty on 
L. Thus the upward fluctuations in L are 
associated with larger weights in the calcula- 
tion of the mean decay length. 

• Radiative events have lower r momenta, 
which tend to result in smaller decay lengths. 
But the lower momenta also tend to yield 
larger opening angles of the daughter tracks 
and therefore smaller uncertainties on L. 
Thus smaller decay lengths are associated 
with larger weights. 

• In some vertex fitting programs, the covari- 
ance matrices describing the errors on the 
reconstructed track parameters are "swum" 
to the location of the fitted vertex and a sec- 
ond fitting iteration is performed. This ap- 
pears to be a reasonable thing to do, but such 
a fitting program assigns larger weights, on 
average, to the r's with larger decay lengths, 
leading to a bias on the average decay length. 

5.3. Detector alignment 

Tracking systems are calibrated and surveyed 
based on tracks reconstructed in the data. This 
procedure is not perfect; after alignment, the av- 
erage impact parameter (d) (which would ideally 
be zero) can vary with 9 and <fi by 10 fim or more. 
The t lifetime measurement is, however, based on 
impact parameters on the order of cr r = 87 /xm. 
How can the experiments claim systematic uncer- 
tainties related to detector alignment as small as 
0.1%? 

A conjecture, put forth by ALEPH fl9|| , pro- 
vides some insight. They theorize that the effects 
of d offsets cancel, to first order, if there are no 
azimuthal holes in the acceptance of the track- 
ing system. To illustrate this point, I present 
the preliminary results of a simple Monte Carlo 
study. Simulated e + e _ — > t + t~ events at y/s = 
91.2 GeV were generated, and a sample of 500 000 
three-prong decays was selected with reasonable 
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Figure 5. Monte Carlo study of impact parame- 
ter offsets, Experiment B (uniform offset plus one 
excursion). The "Input" plot shows the offsets 
applied to the impact parameters d as a func- 
tion of </>. The "Output" plot shows the resulting 
relative bias on the reconstructed decay length, 
B = (L Iec - L tIue ) / (L tru e) as a function of 4>. 



cuts on the momentum and polar angle of the 
daughter tracks. Rather sterile conditions were 
maintained for the experiment: Gaussian track- 
ing errors were generated, with ad — 30 fim, 
= 0.2mrad, and (Sd ■ 5<j)) = O.9cTd0>. A two- 
dimensional vertex fit was performed for each de- 
cay, and no errors on the r production point or 
the r direction were introduced in the calcula- 
tion of the decay length. Impact parameter off- 
sets were then applied as a function of <f>, and the 
fits were repeated, to study the effect on the de- 
cay length bias (L rcc — L trU c)- I tried four different 
d offset configurations, as described below. 

A. No d offsets. When no systematic offsets are 
applied to the impact parameters, the av- 
erage bias is B = (L mc - L tIUC ) / (L tIUC ) = 
(+0.16 ± 0.04)%, reflecting the small posi- 
tive bias due to the correlation of the d and 
(j> errors. The rule of thumb is that the 
relative decay length bias is roughly equal 
to the ratio of the detector-induced correla- 
tion of d and ip to the lifetime-induced cor- 
relation. In this case the detector-induced 
correlation is (8d ■ 5<fi) — 0.0054 /im, while 
the lifetime-induced correlation is roughly 
(d ■ i/j) = (cr T ) (I/7) = 3.4 /im, so the ratio 
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Figure 6. Monte Carlo study of impact param- 
eter offsets, Experiment C (radial shift of silicon 
vertex detector wafers). 



Figure 7. Monte Carlo study of impact parame- 
ter offsets, Experiment D (broken azimuthal ac- 
ceptance) . 



is 0.0054/3.4 = 0.0016, which is comparable 
(!) to the observed offset. 

B. Uniform offset plus one excursion. The "In- 
put" plot in Fig. H shows the applied d offset 
of +20 /im everywhere, plus a triangular ex- 
cursion of amplitude —60 /im in one region. 
The "Output" plot shows that the local bias 
B is essentially the derivative of the input 
function with respect to (f>, with the sharp 
edges smoothed out over an angular scale 
corresponding to the typical opening angle 
of the r decays. In particular, the slope of 
the offset function in the region of the excur- 
sion is ±(60 /im) / (30°) = ±115 /im (convert- 
ing the degrees to radians). This quantity is 
equal to ±6.1% of the mean decay length in 
the xy projection, whereas the maximum ob- 
served bias in the output plot is about ±6%. 
In spite of the large local biases, the global 
average bias is B = (+0.17 ± 0.04)%, i.e., 
unchanged from Experiment A. If the accep- 
tance is unbroken in 4>, and the bias is the 
derivative of the input function, then the av- 
erage bias is proportional to the integral of 
the derivative, which is zero for any input 
offset function. 

C. Radial shift of silicon vertex detector faces. 
Here I consider a one-layer vertex detector 
with nine flat faces at a radius of 6 cm. I 



suppose that the silicon wafers are shifted 
away from the origin by 100 /im with respect 
to their assumed locations. I then make 
the crude approximation that this shift has 
no effect on the reconstructed track direc- 
tions and simply introduces an offset on d 
given by (—100 /im) sin a, where a is the az- 
imuthal angle of incidence of the track on 
the wafer. This scenario corresponds to the 
input function shown in Fig. ^. Again the 
output plot looks like the derivative of the 
input: the local bias has positive S functions 
(again smeared due to the opening angles of 
the r decays) in the regions where the daugh- 
ter tracks do not all pass through the same 
vertex detector face and a uniform negative 
value elsewhere. The bias is locally as large 
as +28%. Nevertheless, the global average 
bias is unchanged: B = (+0.14 ± 0.04)%. 

D. Broken azimuthal acceptance. This experi- 
ment is the same as Experiment C, except 
that 2° azimuthal gaps are introduced be- 
tween adjacent faces of the vertex detector. 
I reject r decays in which one or more of 
the daughter tracks passes through a gap. 
Naturally, this has no effect on the output 
plot (Fig. 0) except near the positive spikes; 
some of the events that had a large posi- 
tive bias in Experiment C are now removed 
from the sample. The resulting negative bias 
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on the global mean decay length is huge: 
B = (-3.73 ± 0.04)%. If the gaps had been 
wider, such that the three daughter tracks 
in selected events always pass through the 
same face of the vertex detector, the bias on 
the mean xy decay length would have been 
simply —100 /im (the radial shift of the faces) 
or -5.3%. 

These experiments show that d offsets due to 
alignment and calibration errors in tracking sys- 
tems have little effect on the measured lifetime in 
the DL method, as long as the azimuthal accep- 
tance is unbroken. In fact it is straightforward 
to prove that the same result holds for the IPD 
method. But these conclusions rely on the as- 
sumption that the weighting of the events in the 
lifetime averaging procedure is independent of <p. 
In a realistic situation where the smearing related 
to the luminous region is significant (not at SLD!) 
and depends on <fi, the (reweighted) integral of the 
derivative of the offset function is, in general, not 
equal to zero. 

It should be mentioned that correlated offsets 
in d and (f> can yield a bias even when the az- 
imuthal acceptance is unbroken. It is straightfor- 
ward to measure the d offsets versus 9 and <fi from 
Bhabha, dimuon, or qq events. Corrections may 
then be applied to the r data. On the other hand, 
offsets in <f> are difficult to measure; the system- 
atic uncertainty on the lifetime should allow for 
a range of possibilities. 

Systematic offsets in impact parameter and di- 
rection may also be present in the rz view. These 
offsets, affecting methods such as DL and 3DIP, 
which do not operate solely in the xy projection, 
are difficult to measure. Moreover, there is no 
bias cancellation rule for such offsets because the 
unbroken acceptance and periodic boundary con- 
ditions do not apply in 6 as they do in <fi. 

Finally, it is interesting to note that the ab- 
solute scale of the lifetime measurements is set 
by the impact parameter scale, which depends on 
the detector dimensions. In experiments with mi- 
crostrip or pixel vertex detectors, it is the pitch of 
those detector elements that counts, not the radii 
of the vertex detector layers. 

6. Treatment of biases 

With some methods, the r lifetime measure- 
ment must be "calibrated" by means of Monte 



Carlo events. For example, the interpretation of 
the impact parameter sum distributions studied 
in the MD/MIPS methods is based on simulated 
events. In other methods (e.g., DL and IPD) 
there is a simple geometric relation between the 
observables and the mean lifetime, so, to first 
order, no calibration is needed. In such cases 
the Monte Carlo is normally used to check for 
"possible" biases in a measurement; a sample of 
events with known generated lifetime is passed 
through the analysis, and the reconstructed mean 
lifetime is compared to the input value. The ex- 
perimenters may then choose one of two valid ap- 
proaches: 

1. If T out put is significantly different from Ti np ut, 
apply the difference as a correction to the 
lifetime measured in the data. 

2. Always apply the difference as a correction 
to the lifetime measured in the data. 



Approach 1 is used by most experiments. It 
turns out that between 1985 and 1993 the dif- 
ference T outp ut — Tinput was not subtracted in 
eight published measurements of r r with the DL 
method. In all eight cases, r ou t pu t was greater 
than Ti n p U t p8j ■ This observation is experimental 
evidence that common biases are present in all 
DL measurements. 

Corrections must be applied for all biases be- 
fore a valid world average can be calculated. I 
would encourage experimenters to apply those 
corrections themselves, but to go far beyond Ap- 
proach 2. We can and should identify and mea- 
sure each individual source of bias in our analy- 
ses by means of Monte Carlo events. Here is an 
example from the IPD method. The bias that re- 
sults from radiative events in which the r + and 
r~ are not back to back in the xy projection can 
be evaluated by calculating A(j> (see Section 2.4 ) 
from Monte Carlo truth information, with and 
without a correction for the acoplanarity of the 
t's, and comparing the fitted lifetimes in the two 
cases. Such a technique yields a measurement 
of this one bias with very good statistical preci- 
sion, and the reliability of the simulation can be 



checked by searching for e + e 



7T T T 



events 



in data and Monte Carlo. 

A series of measurements of this type can be 
devised, using various pieces of information from 
the Monte Carlo truth, in order to evaluate each 
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Figure 8. RPP world average r lifetime value 
versus year Q . 



contribution to r outp ut — Tinput- We can rigor- 
ously evaluate the systematic uncertainty on the 
lifetime measurement only by understanding the 
magnitude of each bias contribution. 

7. History of the r lifetime 

The world average r lifetime values evaluated 
by the Particle Data Group since 1986 are plotted 
in Fig. ||. A fairly steady decline is observed in 
these averages. In my opinion, three factors con- 
tribute to this decline: statistics, unsubtracted 
biases, and "other effects." As evidence for the 
presence of "other effects," the r lifetime average 
in the 1992 Review of Particle Properties Jl4|] had 
X 2 — 2.0 for 10 degrees of freedom, corresponding 
to a confidence level of 0.9962. 

8. Thoughts on the next method 

Here are a few ideas about a possible ultimate 
method for analyzing 1-1 topology events. We 
may be able to squeeze a little more out of the 1-1 
(and perhaps 1-3) topology events if we created a 
method that takes into account all of the available 
information: energies and directions of charged 
and neutral particles, charged particle identifica- 
tion, impact parameters in r<p and rz views, and 
the position and size of the luminous region. 

The method should transform this information 
into two or three variables from which the life- 
time is extracted. It should take into account 
the known r decay dynamics and allow for initial 
and final state radiation. Furthermore, I believe 



we will not be able to make substantial advances 
in precision unless the new method allows the im- 
pact parameter resolution to be fitted from the r 
data itself, as in the 3DIP method. 

9. Outlook and conclusions 

In conclusion, the world average r lepton life- 
time is 

t t = 290.5 ± 1.0 fs, 

and the % 2 of the average looks healthier. The 
LEP experiments dominate the world average. 

It will be a challenge to achieve the next fac- 
tor of 8 improvement in precision on r T . We will 
need to rely on the large statistics of CLEO and 
the b factories, with considerable care to under- 
stand and reduce systematic errors related to the 
fitting procedures, tracking resolution, and back- 
grounds. 
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